Edge excitations and Topological orders in rotating Bose gases 
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The edge excitations and related topological orders of correlated states of a fast rotating Bose 
gas are studied. Using exact diagonalization of small systems, we compute the energies and number 
of edge excitations, as well as the boson occupancy near the edge for various states. The chiral 
Luttinger- liquid theory of Wen is found to be a good description of the edges of the bosonic Laughlin 
and other states identified as members of the principal Jain sequence for bosons. However, we find 
that in a harmonic trap the edge of the state identified as the Moore- Read (Pfaffian) state shows a 
number of anomalies. An experimental way of detecting these correlated states is also discussed. 
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It has been argued Q, H, Q that quantum fluctuations 
can destroy a Bose-Einstein condensate if it rotates very 
fast. Since large amounts of angular momentum can 
be imparted to a cold atomic gas, experimentalists have 
been able to create systems with a large number of vor- 
tices, Ny Thus, the question of what happens 
when iV„ eventually becomes comparable to the num- 
ber of particles, N, has been raised P, 0- In Hi it 
was shown that for v = N/Ny ~ 10, the BEC and the 
Abrikosov vortex lattice are destroyed and replaced by a 
series of "vortex liquid" states, some of which are incom- 
pressible and exhibit good overlap with bosonic versions 
of wavefunctions known from fractional quantum Hall ef- 
fect (FQHE). 

Good overlap is usually a strong indication that some 
of the correlations of a state obtained from exact diago- 
nalization are captured by a model wavefunction. How- 
ever, Wen has emphasized that FQHE states possess 
a new type of quantum order, which he dubbed 'topolog- 
ical order', that provides a better way to classify them. 
In a deep sense, topological order can be regarded as a 
measure of the quantum entanglement existing between 
the particles in a correlated quantum Hall (QH) state . 

In this paper we study the topological order of the 
vortex liquids as reflected in their edge properties 0. 
In the rotating frame, edge excitations are the low-lying 
excitations of the vortex liquid and, contrary to the 
ground states P, |3, O, E|, so far they have received 
little attention. Based on the strong similarities with 
electron FQHE physics, chiral Luttinger liquids 0] and 
similar edge excitations are expected, but an explicit 
demonstration is lacking for a harmonically confined gas 
of bosons under rotation. This is provided here by nu- 
merically diagonalizing the Hamiltonian of small systems. 

The topological order can be studied at the edge of 
fairly small droplets of QH liquid 0, In what 

follows, we shall focus on three types of states that have 
been identified in previous works on small droplets [lllllj|. 
In the high angular momentum end, we study the edge 
properties of the Laughlin state (corresponding to a bulk 



filling fraction 



i). By decreasing the total angu- 



lar momentum, we come across the compact composite- 
fermion states discussed in Ref. Amongst these 

states, we find evidence for the topological order (and re- 
lated edge structures) that correspond to bulk states with 
filling fractions v = ^ and z/ = | of the principal Jain se- 
quence. At even lower angular momentum, we study the 
edge properties of the state identifled0 as the finite-sized 
Moore-Read (or Pfaffian) state^^ (y = 1). We observe 
a number of anomalies that persist up to the largest sizes 
studied here {N = 13). The bulk properties of the states 
at these filling fractions have also been studied in exact 
diagonalizations in edgeless geometries Uliaj which have 
established their interpretation in terms of the Laughlin 
{ly = 1/2), composite fermion {ly — 2/3, 3/4) and Moore- 
Read (zy = 1) states. For the latter, calculations on the 
sphere 0| and on the torus found the correct shift 
and ground state degeneracy (which are also consequence 
of the topological order Q). 

An ultracold gas that rotates rapidly in a cylindrically 
symmetric harmonic trap acquires a pancake shape and 
eventually becomes quasi two-dimensional (2D) when the 
chemical potential fj, < huj^^, where W|| is the axial trap- 
ping frequency Furthermore, the Coriolis force acts 
as an effective Lorentz force, which in a quasi-2D system 
leads to Landau levels (LL's) separated by an energy 
2huj± For ^ < 21iuj±, all atoms lie in the lowest 

Landau level (LLL), and the total single-particle energy 
is proportional (up to a constant) to the angular momen- 
tum. Here we shall be interested in this limit, which has 
been already achieved in the experiments j^]. In the ro- 
tating frame, the Hamiltonian (relative to the zero-point 
energy) is: 



Hlll =T + U2 = Hlu^ - n)L + .9 ^ 5{r, 



(1) 



TiL being the total axial angular momentum and g = 
■\/87r?iw_L(as/^||)^^ the effective coupling for a gas har- 
monically conflned to two dimensions {as < i\\ being the 
scattering length that characterizes the atom-atom inter- 



2 



action in 3D, £|| = yJfi/Muj\\ the axial oscillator length, 
and I — \J1iJMujZ, M being the atom mass). 

Laughlin state: For L = Lq = N{N — 1), the ground 
state of i/LLL is flSini: 



'^>o{zi,...,zn) = {Y[{2 



(2) 



where Zj = {xj+iyj)/£. This wavefunction vanishes when 
any two particles coincide and therefore U2\^o) = 0. 
This property is maintained if it is multiplied by an 
arbitrary symmetric polynomial of the z^. Edge ex- 
citations of the above state are generated Ig by 
elementary symmetric polynomials of the form Sm — 
Sii<i2< - <im Zi^Zi^ ■ ■ -Zi^. Thus, L is increased by m 
units and, according to the excitation energy is 
fi(u± — 0)m. In the rotating frame, the edge excitations 
are the lowest energy excitations of a rotating Bose-gas in 
the Laughlin state 8r| . Their degeneracy for L = Lq + m 
is given by the number of distinct ways m can be written 
as a sum of smaller non-negative integers {i.e. partitions 
of TO, p(rn)). For example, for m = 4 there are five de- 
generate states: (s4, S3S1, (s2)'^, S2(si)^, (si)**) X $0- The 
properties of these wavefunctions are captured by an ef- 
fective field theory 0,0,0, which treats these excitations 
as a non-interacting phonon system with Hamiltonian 
-ffcdgo = h{uj± - Q)L, where L = Lq + X]m>o ™ ^Lbm, 
and [bm,bl^] = 6mn, commuting otherwise. Hence, the 
number of edge states (NoS) for a given to can be ob- 
tained from a generating (or partition) function: Z{q) = 

Tr g^-^" = n™>o(l - 9™)"' = Em>oPi^)l"'- In ta- 
bleQI we compare the theoretical NoS (A'' —^ 00) with the 
numerical results, finding excellent agreement for m < N. 
The deviations are due to finite-size effects and can be 
accounted by a 'truncated' function Z^{q), where the 
product is restricted to < m < A^. 

Another prediction of the theory 0, 0] the form of 
the groundstate boson occupancy, n(l), just below the 
highest occupied orbital, l^ax — 2(A^ — 1). The follow- 
ing ratios of n(/)/n(Zmax) are predicted: 1 : 2 : 3 : 4 for 
I = Ims^x, ■ ■ • , ^max — 3. For N = 7 bosons in the Laughlin 
state, we find 1.0 : 2.0 : 2.9 : 3.5, which is in good agree- 
ment with the theory given that its validity is limited to 
^max — I ^ ~ 3 13 • Similar agreement was found for 
the fermion Laughlin state in Ref. 0- 

Principal Jain sequence: The incompressibility of 
the Laughlin droplet can be regarded as a consequence 
of a statistical transmutation: at large Ny, a boson 
binds one vortex and becomes a composite object (called 
composite fermion, CF) [lll | that behaves as a spinless 
fermion. In a Laughlin droplet, the CF's fill up the 
N lowest angular momentum orbitals in the LLL. Com- 
pressing the droplet and therefore decreasing L requires 
promoting CF's to higher LL's and costs a finite amount 
of energy. Thus, states with lower angular momentum 
will contain CF's in higher CF LL's, and this leads to 



TABLE I: Number of edge states (NoS) vs. excitation angular 
momentum m. In the data for the Moore-Read (Pfaffian) 
state A = for a pure three-body interaction and A = 1 for a 
pure two-body interaction. 
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Spectrum of the state L = 30 with N = 7. The 
: 0.091 5/^^. Below the gap, we observe two phonon 



FIG. 1: 

gap A = 

branches along with their multiphonon excitations. 



the ansatz: 

$(zi,...,Z7v) 



V[lliz,-z,)'i>^g^^\r^,...,r^)]. (3) 



and is a 



where V projects onto the LLL 11] 
Slater determinant with {Ni}^^^ CF's filling the lowest 
angular momentum orbitals of p > 1 CF LL's and L = 
Lo = N{N - l)/2 + NiiN, - {2t + l))/2 (No = 

N and Ni = for i > in the Laughlin state). In 
what follows, we focus on the CF states {4, 2} and (5, 2}, 
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TABLE II: Interaction energies (in units of g/i^) of the edge 
excitations of the {5,2} state (A'^ = 7, L = 30, see Fig. 0. 
The predictions for multiphonon states are given in brackets. 
Deviations are due to non- linear terms not included in Eq. Q . 
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which have been shown 1 1] to have good overlap with the 
exact states at Lq = 20 {N = 6) and Lq = 30 (iV = 7), 
respectively. For a relatively large range of L — Lq > 
(3 for {4, 2} and 4 for {5, 2}) the lowest energy state is 
a centre-of-mass excitation of the state at L = La, and 
the interaction energy is unchanged. Edge excitations 
are those states with energy lower than the bulk gap (A, 
cf. Fig.Ul) [3. According to [3, ISj , these states exhibit 
two branches of edge phonons described by: 



Hc.dnc. — Ef) 



m>0 



(4) 



m>0 ' 



with [dm, dl] 



and L - 

= 6nm, Commuting otherwise. Note that the 
edge excitations are not degenerate in energy. However, 
one can still compute the NoS from Z{q) — Tr q^^^° — 
T\m>oi^ — q™)~^. In table^we compare the numerical 
results to the theoretical predictions for the NoS, finding 
perfect agreement for m = 1,2. For higher m the results 
are affected by finite-size effects. However, in the state 
{5, 2} the NoS for m = 3 is quite close to its TV — > oo 
value, and therefore we concentrate on this state for fur- 
ther analysis. We next try to reproduce the energies of 
the sixteen edge states of {5, 2} from m = 1 to 3 using 



,(6) 



L2 s-nd (^m=i 2 the only fitting parameters {tOi 
has only a kinetic energy contribution by Kohn's theo- 
rem 1^). The interaction part of hujf^ can be extracted 
from the m — 1 data: Tiwj'^^ = 0.015 For m = 2, 

there are the five following states: 2-'^/^ (&I)^|0), ^lO)' 
d\bl\0), 2-1/2 (4)2|o), 4|0), whose energies can be ob- 
tained from Eq. I^J. For m — 3 the states and energies 
can be written down in a similar fashion. The compari- 
son with numerics for to = 2, 3 is given in table UTI Thus, 
these states are good representatives of the topological 
order of the filling v = \ from the principal Jain sequence 
for bosons 0, Q . 

Finally, we also find evidence (to be reported else- 
where |22j ) for the edge structures (three phonon 
branches) corresponding to v = |. However, the NoS 
for m > 1 is affected by finite-size effects for N = 6, 7. 



Moore-Read (MR) or PfafHan state: Ground 
states of CF's with lower angular momentum are ob- 
tained by placing CF's in higher effective LL's. Even- 
tually, when the number of occupied levels p — > -|-cxi, the 
CF's would not feel any effective Coriolis force and the 
resulting state should be compressible. However, in such 
a state the CF's can pair and condense into a BCS state, 
which would render the state incompressible again |l5l |. 
Since CF's are spinless, the pairing takes place in p-wave 
and the BCS wavefunction is a Pfaffian [ia, ll^ , which 
must be multiplied by rii<j(-^» ~ ^j) yield a bosonic 
wavefunction: 



MR 



i<j 



Zj)A 



1 



Zl - Z2 



zn-1 — zn 



(5) 



where A stands for anti-symmetrization (see e.g. Refs. P, 
ll^ ) of the bracketed product, which is the Moore- Read 
(MR) wavefunction j 1 6| . Soon after the MR state was 
introduced, it was pointed out jl^j that it is a zero 
interaction-energy eigenstate of a three-body potential, 
U3 = 93Y.i<j<k^i^i ~ - Tfe)- Thus, it is con- 

venient to work with a modified Hamiltonian: H^j^^ — 
T + AC/2 + (1 — A)C/3, so that for A = the ground state 
at angular momentum ~ N{N — 2)/2 (even N) or 
^{N- 1)2/2 (odd N) is exactly the MR state. In- 
deed, this modification is not entirely artificial as one of 
us has recently shown that MR can become exact near a 
Feshbach resonance |l^ . 

For the exact MR state (i.e. A = 0) besides the polyno- 
mials Sm. introduced above. Wen and Milovanovic and 
Read [ig found a branch of fermionic edge excitations 
which are generated by replacing the Pfaffian in $mr by 

A [z"^' ■■■zl" {ZF+I - ZF+2)~^ ■ ■ ■ {ZN-I - zn)~^] , 

where ni,...,np are non-negative integers. 
Thus the angular momentum is increased by 
L — = ('^fc + \)- This spectrum, together 

with the phonon branch related to Sm is described 
by L = L^R + Y,^^^^ 5t^5^ + (m - i) c|^_ 1 c,„_ 1] , 

where hm, are the phonon operators, and the fermions 
{c„j_ 1 , ~ 5mm anti-commuting otherwise. How- 

ever, due to the paired nature of the state, even and 
odd A'' are different. For instance, to compute the NoS, 
one must define Id Z'=™'^(g) = iTr (1 + {~lY)q^~^° 



and Z°'^'^ 



[q] = iTr (1 - (-1) 



' )q since the parity 

(-1)^\ with F = Em>oc]„_iCm_i , is a good quantum 
number 6]. The numerical results are compared with 
the NoS tor N ^ 00 in table Perfect agreement is 
found for to < [Af/2]. Higher values of to are shown to 
illustrate the effects of finite size; the observed differ- 
ences from the A^ — > 00 values can be also accounted 
for by the theory Furthermore, using the effective 

field theory we have also obtained the behavior of 
n{l) near the edge. For A'' even the predicted ratios of 



n{l)ln{lmax), where 1^ 



A^ - 2, are 1 : 2 : 3 : 



For 



4 



2.0 




Orbital angular momentum I Orbital angular momentum 1 



FIG. 2: Orbital occupancy for the Moore-Read (or Pfafiian) 
state (3-body) and the eigenstate of the 2-body interaction 
identified as the Pfaffian. 



= 8, we numerically find 1.0 : 2.2 : 3.4. However, for 

N odd the ratios of n{l < Zmax)/"(^max) ('max = - 1) 

behave differently 1:1:2:3. For iV = 9 we 

numerically find 1.0 : 1.0 : 2.5 : 4.1, and the scaling 
observed from smaller systems shows a trend of conver- 
gence to the predicted ratios. Although n{l) had been 
analyzed in |6| for a fermion MR state, the different 
behavior of Ti{l) for odd N had not been described. 

As soon as the two-body interaction is turned on, i.e. 
already for small A > 0, we observe that the plateau 
at L > L^^ is lost for small L - L^^. For pure two- 
body interaction (A = 1), the states at + 1 for 
TV = 5, 9, 11, 13 and at L^^ + 2 iov N = 6, 7, 8, 10, 12 
have lower interaction energy than the ground state at 
j^MK^ Thus, U2 strongly perturbs the plateau of the ex- 
act MR state by favoring states with angular momenta of 
compact CF states. If, regardless of this fact, one counts 
the number of states at L > with interaction ener- 
gies less than the first excited state at Lq^^, the results 
do not seem to converge and disagree with the NoS ex- 
pected for the MR state (see table Qlfor the observed NoS 
at A = 1 for N = 12, 13). We also observe a rapid deteri- 
oration of the overlap of the state at with the exact 
MR state: from 0.91 {N = 5) and 0.90 {N = 6) to 0.68 
(TV = 8) and 0.74 (N = 9). These discrepancies might 
be due to a very slow convergence as N grows towards a 
well-defined bulk MR state. However, we note that for 
very large N, local density arguments |21| show that the 
edge of the MR state will reconstruct. The anomalies in 
the edge of the small systems observed here are striking 
in view of the good behavior exhibited by the Jain states, 
which are also approximate wavefunctions. 

Experimental consequences: It is possible to excite 
the surface modes by inducing a small time-dependent 
deformation of the harmonic trap. Within linear re- 
sponse, the energy injected by an m-polar deformation is 
proportional to the oscillator strength /„j = |(Lo + 
m, a|Om|Lo)P, (Om = Y^iLi zD, where a runs over aU 



edge excitations. The theory predicts |1| /m = mvE?'^, 
where R = yjNjv is the droplet radius. For the dipole, 
/i — N, as required by Kohn's theorem 0, and con- 
firmed by our numerics. For the quadrupole (m = 2) 
/2 = 2N'^/iy. We have numerically tested the accuracy 
of this formula for the m = 2 modes of the states in 
table m (except the MR state at A — 1), finding that al- 
ready for iV = 5, 6 the deviations are no larger than 7%. 
This result suggests an experimental way of estimating 
the filling fraction, provided /2 and N (or /i) can be 
measured: = /2/2iV^. More details will be provided 
elsewhere |22l |. 
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